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Abstract 

We show that recent data from HERA on the proton structure function F2 at small x and large Q 2 
provide a direct confirmation of the double asymptotic scaling prediction of perturbative QCD. A 
linear rise of In F2 with the scaling variable o is observed throughout the kinematic region probed 
at HERA, and the measured slope is in excellent agreement with the QCD prediction. This 
provides a direct determination of the leading coefficient of the beta function. At large values of 
the scaling variable p the data display a small but statistically significant scaling violation. 
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Perturbative QCD predicts that at sufficiently large t = lnQ 2 /A 2 and small x the 
nucleon structure function F 2 should exhibit double scaling in the two variables 



provided only that the nonperturbative input to the perturbative evolution is sufficiently 
soft. We have shown QTJ] that this prediction is indeed confirmed by the first measurements 
of F% performed at HERA0||]. In fact, it turns out that not only most of the HERA data, 
but even some of the older data from the NMC [|J, lie well inside the asymptotic regime, 
suggesting that the starting scale to = hiQoV-^ 2 ^ or ^ ne perturbative evolution should be 
little more than Qq ~ 1 GeV 2 . A significantly enlarged set of measurements of F% has now 
become available ||[§, which makes it possible to test double scaling more quantitatively 
Specifically, the slope of the linear rise of lni 7 ^ in the scaling variable a can be reliably 
measured, and turns out to be in excellent agreement with the QCD prediction, thus giving 
a direct empirical determination of the leading coefficient (3q of the QCD beta-function. 
We also find that there is now evidence for scaling violation at large p. 

Double asymptotic scaling follows from a computation of the asymptotic form 
of the structure function F^x^t) at small x based on the use of the operator product 
expansion and renormalization group at leading perturbative order. It thus relies only on 
the assumption that any increase in F% (x; t) at small x is generated by perturbative QCD 
evolution, rather than being due to some other (nonperturbative) mechanism manifested 
by an increase in the starting distribution F^x^to). The resulting asymptotic behaviour 
takes the form 



where 7 = 2a/A^ c //3 , (3 = ^-N c - |n/, 6 = (1 + - ttn;), and the unknown 



function /, which depends on the details of the starting distribution, tends to one for 
sufficiently small values of its argument. is an a priori undetermined normalization 
factor. 

In |jj we derived @ by noting that at small- a; the one loop QCD evolution equations 
reduce to wave equations, which propagate the parton distribution functions from their 
boundary values at t = to an d x = xq to larger values of t and smaller values of x. Since 
the propagation is unstable, away from the boundaries an exponential increase with a of 
the form (|]) inevitably arises, provided only that the small-x behaviour of the starting 
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distributions at to is sufficiently soft (which in practice means that if f s (x;t) is a singlet 
parton distribution function, x 1+x f s (x; to) — > as x — > for any A ^ 0.2). The behaviour 
(H) is thus a rather clean prediction of perturbative QCD, in so far as it is independent of the 
details of the (soft) nonperturbative parton distributions which are input at to, provided 
that at small x these conform to expectations based on Regge theory. The asymptotic 
behaviour can be shown [0J to set in rather rapidly as a increases in a region not too close 
to the boundaries, i.e. when p is neither too large nor too small. 

In order to compare the data for F% with the prediction (0) we rescale the measured 
values of by a factor 



to remove the part of the leading subasymptotic behaviour which can be calculated in a 
model independent way.0 Then In [R' F F2\ is predicted to rise linearly with cr, independently 
of p (when p is large), with slope 



if rif — 4 as in the HERA kinematic range. The model-dependent subasymptotic behaviour 
due to the function / can be eliminated by cutting all points with subasymptotically small 
p; the scaling analysis of Ref.[l] (see fig. 2 below) suggests that we place the cut at p 2 = 2. 

All the available experimental data0,|[|| for which pass this cut are plotted in 
fig. [l]. The predicted linear rise in a is spectacularly confirmed, providing clear evidence 
that in the region a 2 > 1, p 2 > 2 the asymptotic behaviour (§) has set in. Indeed, the 
scaling actually sets in rather precociously: even the NMC data down to a ~ 0.7 seem to 
be rising linearly, with possibly an indication of a systematic normalization mismatch of 
around 10% between the NMC and the HERA determinations of F^. 

Fitting a straight line to all 80 HERA points in the plot yields a x 2 of 66, and a gradient 
27ex P = 2.37 ± 0.16, in perfect agreement with the QCD prediction eq.(|]). Turning this 
into a measurement of the leading coefficient of the beta-function gives (with N c = 3) 
(3q = 8.6±1.1 (to be compared with 25/3 for nf = 4). This is a direct, model independent, 
and highly nontrivial test of the perturbative dynamics of asymptotically free nonabelian 
gauge theory. 

1 The constant rescaling factor R may of course be chosen arbitrarily; here we choose R = 8.1, 
so that the normalization of the figures is the same as in O. 



R' F (a, p) = iZexp (5(cr/p) + | In a + ln(p/ 7 )), 
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N s 


N h 


x 2 


a) 


0.341 ±0.005 





96 


b) 





0.156 ±0.002 


878 


a)+b) 


0.319 ±0.012 


0.012 ±0.002 


91 



Table: The fitted normalizations N s and Nh and the associated x 2 s (103 data points). The 
different cases considered are a) soft pomeron b) hard pomeron, and the linear combination 
a) + b). 

We next consider scaling violations, both in the subasymptotic region of small a and 
small p, and in the post-asymptotic region of large p. This is best done by rescaling 
by a factor 

R F (a, p) = Rexp ( - 2 7 cr + 5(a/p) + \ In a ± ln(p/7)) (5) 

to remove all the leading behaviour in (Q). The rescaled structure function should thus 
scale in both a and p when both are sufficiently large to lie in the asymptotic region: 
RfF^ = N ± 0(1/ a) + 0(1/ p). This double asymptotic scaling behaviour is tested in the 
two scaling plots fig. [|, where we also display the predictions obtained by applying the 
leading small-x form of the evolution equations to a typical soft starting gluon distribution. 
Specifically, fig. |2|a) shows that the scaling in a sets in very rapidly, as all the points on 
the plot lie in the asymptotic regime; fig. 0b) shows that the scaling in p only sets in for 
However even if p is as low as p ~ \ the subasymptotic corrections due to /(7/p) 
seem fairly well accounted for by the scaling violation displayed by the curves of fig. [| 

More interestingly, at large p there now appears to be a statistically significant rise 
above the scaling prediction. To test the significance of this rise, we fitted to the data a 
linear combination of the behaviour discussed above and displayed by the curves of fig. [| 
and a "hard pomeron" behaviour, which violates scaling by rising with p (see ref.jl] for 
a more detailed discussion). Including in the fit the 103 HERA points with both a 2 and 
p 2 greater than one half gives the results displayed in the table. The data seem to prefer 
a 4 ± 1% admixture of the hard pomeron solution. One should be very cautious about 
taking this as evidence for the hard pomeron per se, however, since higher loop corrections 
should give a similar rise||. It should be possible to settle this issue decisively when a 
more detailed set of data and more accurate theoretical calculations become available. 

Acknowledgement: We are very grateful to Dick Roberts for communicating to us 
the data of refs [|5||J] . 
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Figure Captions 



Fig. 1. Values of R'pF^ plotted against a: diamonds are ZEUS data ||, squares HI data 
0, and crosses are NMC data. The best fit straight line is also shown. 

Fig. 2. RpF^ plotted against a) a and b) p. Included in the plots are all the HERA data 
with p > 1.2, a > 0.7, respectively. The curves show the prediction obtained 
evolving a typical soft starting gluon distribution: a) dot-dash curve, p = 1.4; 
solid curve, p = 2.2; dotted curve, p = 3.2. b) dot-dash curve, a = 1.1; solid 
curve, a = 1.8; dotted curve, a = 2.1. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406385v3 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406385v3 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406385v3 
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